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This paper proposes a feedback linearization control of doubly fed induction 
generator based wind energy systems for improving decoupled control of the 
active and reactive powers stator. In order to enhance dynamic performance 
of the controller studied, the adopted control is reinforced by a fuzzy logic 
controller. This approach is designed without any model of rotor flux 


estimation. The difficulty of measuring of rotor flux is overcome by using 





high gain observer. The stability of the nonlinear observer is proved by the 
Lyapunov theory. Numerical simulations using MATLAB-SIMULINK 
shown clearly the robustness of the proposed control, particularly to the 
disturbance rejection and parametric variations compared with the 
conventional method. 


Keywords: 


DFIG 

Feedback linearization control 
Fuzzy logic 

High gain observer 


Lyapunov stability This is an open access article under the CC BY-SA license. 





Corresponding Author: 


Kada Boureguig, 

Oran University of Science and Technology - Mohamed Boudiaf, 
USTO-MB, BP 1505 El M'Naouer, 31000 Oran, Algeria. 

Email: kada.boureguig@univ-usto.dz 








1. INTRODUCTION 

The wind energy is a pollution-free and effective source. Therefore, a wind power generation system 
becomes one of the potential sources of alternative energy for the future [1]. Energy consumption over the 
last century has increased significantly due to the great industrialization. Recently, particular interest has 
been given to generating electricity from renewable energy sources. Of all renewable sources, wind energy 
holds the largest market share and is expected to maintain rapid growth in the coming years [2]. Wind energy 
systems have received considerable attention over the past decade as one of the most promising renewable 
energy sources due to negative environmental influences and the high cost of conventional energy sources. In 
this context, several countries have turned to explore the wind energy sector, leading researchers to conduct 
research to improve the efficiency and power of electromechanical conversion and quality of providing 
energy [3].To meet energy needs, it is imperative to find adjusted and flexible solutions by reducing energy 
consumption or increasing energy production by adding power plants or improving the efficiency of existing 
installations. In addition, the dynamic improvement of the performance of renewable energy systems, whose 
non-linear characteristics are particularly important, especially with the rapid growth of their use. Therefore, 
the control of the Wind Power Conversion System (WECS) based on doubly feed induction generator 
(DFIG) with intermittent input wind speed is particular interest in the energy and control communities. 
DFIGs are potential candidates for high power wind systems because they can generate reactive current and 
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produce constant frequency power for variable speed operation. However, the main disadvantage of the 
DFIG is its brushes and slip rings structure, which involves permanent maintenance and reduces the life time 
of the machine [4]. This type of generators has been widely used for wind systems [5]. The control and 
operation of these systems has been the subject of many research projects in recent years [6]. In this 
perspective, various studies have been proposed in [7-9] to control wind turbine systems based DFIG from 
the classical configurations by application of Field oriented control (FOC) strategy. However, the WT-DFIG 
is highly nonlinear system, with strong couplings between the different variables of the systems. 

In this context, many nonlinear control methods have been developed. Penghan Li et al [10] 
proposes a non linear controller based on state feedback linearization strategy to reduce sub-synchronous 
control interaction inseries-compensated doubly fed induction generator (DFIG)-based wind power plants. 

In [11] the authors have used a robust nonlinear feedback linearization controller based sliding mode 
control to relieve sub-synchronous control interaction in doubly-fed induction generator based wind farms 
connected to series-compensated transmission lines. The two aforementioned works show good 
performances. However, the authors did not use a nonlinear model and DFIG Wind Turbine control scheme 
is based on vector control. In addition, the states of the DFIG are supposed to be measured. 

Djillali et al [12] have used Neural Input-Output Feedback Linearization Control. The neural 
controller is based on a Recurrent High Order Neural Network, trained with an Extended Kalman Filter. This 
last method uses a simple PI controller to define the control law defined by a relationship linking the new 
internal inputs to the physical inputs. The same strategy control has been applied in [13] based on the linear 
quadratic regulator (LQR). Due to the limitations presented by the two linear controllers the PI and the LQR 
which have a low robustness to parameter variations as well as to unbalanced grid voltage which have direct 
effects on the dynamic performance of the system, and poses serious problems, such as oscillations of the 
stator power and the generator torque, which are detrimental to the mechanical system and the electrical 
network [14]. 

In this paper the Feedback linearization technique is combined with Fuzzy logic to form Fuzzy- 
Feedback Linearization Controller applied to a non-linear model of DFIG to improve the performance of the 
system such as the response time , robustness against parameter variations and the sensitivity to perturbations 
(unbalanced grid voltage) .This new method is augmented by High Gain Observer (HGO) mainly used to 
estimate generator rotor flux, based on the measurement of rotor currents ,stator voltages and the mechanical 
speed. The effectiveness of the proposed controller is compared to the conventional Feedback linearization 
control by simulation results in Matlab Simulink. 


2. MODELING OF STUDIED SYSTEM 

As shown in Figure 1 the system is composed of two parts; the first is the conversion of the kinetic 
energy of the wind into mechanical energy via a turbine and the second is the conversion of the mechanical 
energy at the level of the turbine shaft into electrical energy via a double-feed generator. The stator is 
connected directly to the grid and its rotor also via a static converter which allows delivering the necessary 
control voltages of the stator powers. 
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Figurel. System under study 
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2.1. Turbine model 
The aerodynamic power appearing at the rotor of the turbine is then written [15]. 


1 
P; = +pCp(A, B)SV? (1) 
Where: p is the air density and S the swept surface area of the turbine (xR?) 
V is the wind speed (m/s), Cp(A,B) is the power coefficient of the turbine, A is the tip speed ratio 
and $ is the pitch angle. The tip speed ratio is: 


= at (2) 


where: R is the radius of the turbine (m) and œ+ is the speed turbine (rad/s).Figure 2 shows the curve of the 
power coefficient versus A for a constant value of the pitch angle f. 
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Figure 2.Typical curve of power coefficient. 


2.2. Dynamique model of DFIG 
The mathematical models of three DFIG phases in the Park frame are written as follows [16]: 


: dos 
Vas = Rsids + at Os qs 


dọ 
= as 
Vas = Rsigs + = + OsQy, 





es (3) 
Var = Rriar + dt (®, = ©) gr 
> War 
Var = Rrigr + G+ (@, — ©)94, 
Pas = Lsias + Miar 
Pqs 7 Lsigs + Migr (4) 


Gar = Lrigr + Migs 
C= Lrigr + Migs 


Where R,, R,are the rotor and stator resistances, respectively; Ls, Ly, M are the rotor, stator and 
mutual inductances, respectively; igs, igs-rare the stator-rotor current components; Vas-r , Vgs—rare the 
components of the stator-rotor voltage; M4... Pqs—rare the stator-rotor flux components; w , œs are the 
rotating and stator pulsations, respectively. 

The model of DFIG according to the rotor components is represented by the following equations [15]. 





digr A . 

at tar + Osigr + 429 gr 430g, — A4Vds + a3Var 
di 

qr _ . 2 

a —Oslgr—Ayzigr + a2Pqr t430P4r — a4Vgs + a3Vqr 





5) 
dog x ( 
E = —Rrlar + OsP gr — OP gr + Var 
do 
Top: — 
a Rrigr OsQ 4, + OP a, + Var 
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where 

1 1 1 1 1—0 1-M2 

aim Gtm a As ae 1- LsLr »b=R,, 
2 P L L 

G=>,0G=+, eta 

1 j 2 j 59S Rs or Ry 

o is the dispersion coefficient 

The electromechanical dynamic equation is then given by 

do _ P? p P 

at je (Qalar g Parlar) F 7 Cvis + Cg) (6) 


where P is the number of pole pairs; J is the inertia of the shaft , Cg is the torque on the generator .All 
frictions on this shaft are included in Cyis. 
We put 


idr iq ® 1 Qa O) = (Zp X2 X3 X4, X5) 
dr’ Yqr 


The system (6) is then written in the form: 


x = f(x) + g&)u (7) 
Where 
d 
= = AE) + a3Vdr 
d 
<2 = fy (X) + a3Vqr 
d 
f(x) =) S38 BE) + var (8) 
d 
St = fax) + Var 
dxs _ 
ETA fs) P 
= T. _[a 0 1100 
u= [var Var] , g(x) = 0 az 0 1 0 (9) 
And 
fı (x) = —a1X1 + OsX2 + a2X3 — a3X5X4 — a4Vds 
f2 (X) = —0sX1 — a1X2 + a2X4 + a3X5X3 — a4Vqs 
f3(x) = —bx, + @,X4 — X5X4, f4 (X) = —bx2 — @5X3 + X5X3 


fs &) = Cy (%4X1 — X3X2) + C2 (Ca — Cyis) 


3. FEEDBACK LINEARIZATION CONTROL 

To develop nonlinear control of the active and reactive powers of stator, the feedback linearization 
strategy is proposed. This technique consists to transform nonlinear systems into linear ones, so that linear 
control techniques can be applied. This technique is possible through change of variables and by choosing a 
suitable control input [17].According to the model of the DFIG developed above, and recalling that the 
reference is chosen so that its component (d) coincides with the stator voltage vector, this system has as input 
variables the voltage applied to the rotor (VarVqr) and as output variables the active and reactive power at the 


stator (P;,Q,) defined by: 


iG = Vasigs + Vasias 
Qs = Vasids — Vasigs 


(10) 


The stator powers’ control law is computed according to the rotor current measurement and 
estimated rotor flux. The latter comes from the proposed high gain observer. Substituting ig, and ig, in (10) 
by their counterparts extracted from the two last equations of (4), one has: 


gr—Lrlqr Og-—Lrla. 
Ps = Vgs = T vs ŒS) 


M M 
11) 
ar —brla 9gr—Lrlqr ( 
Qs = Vąs( ar 7 =) — Vas (~ M ) 
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Arranging (11) 
— Par Lrlar y ar ar 
Ps = -r Vas Tm Vas Hm Yds 
Pa _ Lrlar ® qr rg r 
Qs = aY qs M Vąqs 7 m Vds T Vds 
Differentiating (12) until an input appears 
5 — X4 LyX2 X3 LrX1 
Ps = ria M |, Vas $ mM Vds — M Vds 
LyX X4 “me 
Q; = TM Vas 7 M Vgs — m Vds +t Vds 
Write the last equation as follows 
5 (f3—Lrf1) (f4—-Lrf2) (1-a3Lr) (1-a3Lr) 
P, = m Vas + m Vas + VasVar + VqgsVar 
; me f4) (fg-Lyf1) (1-a3Lr) (a3Ly—1) 
Qs = Vas + M Vas + M VqsVdr + Vas Var 
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(12) 


(13) 


(14) 


It is desired to regulate the output quantities P,and Q,to their respective reference value Poref and Qsref 


For this purpose defining the adjustment errors 


{ e1 = Pref — Ps 
e2 = Qsref — Qs 


Defining the input of the DFIG system 
T 
u= [u, ug]? = [Vor Var] 


Rewriting (14) in the matrix form 








P Ghr y+ (f4-Lrf2) ay, (i=a3Lr) awy a- asly v 
s|— M s s “E 
H 7 Urf) ys + (f3-Lrf1) ah abe“ a- slay 
M S 


Rewriting the new inputs V,and Vin the form (18) 


F =V; 
Q; = V2 
From (19) we can write 
1 
[>] a v] 
Qs 0 $ V2 
S 


From (17) and (18) the control law is given as 


[a] = 07-a +] 


where 
f3—Lrf. f4—Lrf; 1-a3L 1-a3L 
oo Pyy + By C 3 Dys C 3 r) 
= and E(x) = 
A(x) Mrt) y y (ao lef) d ECx) Gsl- Gase) 
M ds M qs M ds M 
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The reference active power Ps is generated by MPPT and the reactive power Qs is defined by the 
grid to support the network voltage .To follow the trajectory of Peer and Qšref , we use a PI controller 
imposed to the linearized system [17]. The new input v is given by 


Paref = kp161 = Kia Í e,dt 


ae hee (21) 
Qsref ~~ Kp2€2 = kiz f ezdt 


V2 





4. HIGH GAIN OBSERVER FOR FLUX ESTIMATION 
We intend to construct such an observer, based on the measurement of the rotor currents, speed and 
voltages [15] the principle is shown in Figure 3. 


[var; Var] 







liari tar] 
[Vas; Vas] So 


[Gar i Şar] 


Figure 3. Rotor flux observation strategy. 


The estimated flux components (@ ar Pgr) are used in the computation of the FBL law.We propose 
to implement a rotor flux observer in order to study its properties. 
From the model (7) and assuming the constant velocity (Q = 0), with Q = > we can write: 


X = A(Q)X + Bv (22) 
where 
X= lira irq Gar Par (23) 


A(Q), B and v are given by 


v= [vsa Vsq VraVrql x (24) 
~a, Os az —pa3Q —a, 0 a; 0 
_ | ®s ay pa3Q a2 p—|0 a, 0 a3 

AD) =| R, 0 0 &—polr’ lo. 0 1 0 25) 
0 =R, —(@, — pQ) 0 o O 0o 1 


Thus, at constant Q, the model is linear, which is a particular case of the form of injection of the 
output and of the output derivative. 

In this part, we are interested in the work presented in [18, 19] which deal with the synthesis of 
observers with high gain for locally observable systems. 

Then it is possible to make out the following change of variables: 


hy 
m = 26 
a oo DS l ) 
z =I 27. 
P = —(aT + 03)! + +(a2 T + aspQ3)o a 


For these changes, model (7) takes the following form: 
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Zi = Z —agV, + a3Vr 
Z) = —(a2T — ©) (Zz — a4Vs + a3V,) + 


os 28 
(a2 T + agpQS)[-R-2 C8) 

—(@s — pQ)S(az T + agpQ3)“*[z. — (aT + @s3)z1] + vr] 
We put 

_ [Z2]  _ [ira sjt] o0 %9 O 71 0.a T0 1). 
E e E a a a ol: 

yis Vr, 9, Z1) 
Ys Vr Q, Z) Ra Ve Q, Z1) Z2) (29) 
i Vp Q, Z1) = —ą4Vs + a3Vr (30) 
Wels, Vr Q, Zi» Z2) = Z2 


That transforms the nonlinear system (28) into a local system of pyramidal coordinates 


f = AZ + W(V5, Ve Z) (31) 


y= Cz 


With the output vector C = [F, 0] 
Then the following system 


2 = A? + yY z) — SZtCT(C2 — y) (32) 
Is exponential observer of the system with Sg is the matrix defined by 


e's 97° 
Sa ee 20-3 G2) 


With 0 > 0 Is the unique solution of the following Lyapunov algebraic equation: 
OS + ATS, + SoA = CTC (34) 


4.1. Theorem 
The function y is globally Lipchitzian with respect to z uniformly with respect to v, and v, 


IVC, Vs, Vr) — WG Vs vel] < £112 — zll (35) 


4.2. Proof of stability anaylsis and observer convergence 
Consider the error 


e=2-z (36) 
Its dynamics is given by 

è = (A — SẸ +CTC)e + y(Z, Vs, Vr) — YÈZ, Vo Vr) (37) 
Let’s consider the following Lyapunov function candidate 

V(e) = e'Sye (38) 
Its derivative is 


V(e) = eTSge + eTSpé (39) 
= [eTa — S§tCTC)T + (vE, Vg Vr) — WG, vov:))] Spe + 
eTSo[(A — Sg*C Ce + yw, vs vr) — YZ Vs Vr)] 
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= eT [A™Sy — 2CTC + SpA] + 2e7So[w(Z, Vs, Vr) — WZ Vs Vy] 
V(e) = —eT (0S + CTC)e + 2eTSp[W(2, Vs Vr) — YC Vs, Vr) (40) 
Using inequality CTC > 0 and theorem (1) we can increase V(e) as following: 


V(e) < —eTASye + 2heT Spe 
< —(0 — 24)V(e) (41) 


This guarantees the exponential stability of the observer for 0 > 2% this concludes the proof. 


4.3. Observer in the initial coordinates 
The observation £ of the state x for the model (7) is obtained by: 


a(z) (42) 


That z = O(x) imply = = ae z 


Is a method of synthesizing the observer expressed in x 





= f(x) + g(x, v) - op S(0)-*CT (CR — y) (43) 


Or in the panoramic form: 


=> 
a 
= 


—aylar + Oslqr t 29 4,—43PQP,, 


* : 7 pe Ji —a, 0 a, O] [Yds 
lgr = —Oshar—Arigr ta3pQG,,. + a29,, 0 -a, 0 a3]]VYaqs E 
Parl =R rfar — (@, + PQVG,, 0 : 1 | vee 
Far —Ryigr + (—@, + pQ)@,, 0 0 oF 
1 0 0 o0 “1720 0 
0 1 0 0 02 20) Par — a 
ay Os a2 —a3pQ lke i “| i 
“Os ~a, —agpQ az @? 


The DFIG control scheme using conventional feedback linearization (CFBL) associated with HGO 
is shown in Figure4. 
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Figure.4.DFIG control scheme using conventional feedback linearization control (CFBL) 
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5. FUZZY-FEEDBACK LINEARIZATION CONTROL 

Fuzzy logic control has been widely used in recent years due to its simplicity of implementation and 
its ability to control nonlinear systems, which gives better performance under parameter variations and 
voltage disturbances [20]. 

Based on the study of fuzzy logic control described in [21], [6] we will proceed to its application to 
DFIG based on the equation (19) or we consider our system as a simple integrator which facilitates the 
synthesis of this control where we will have two regulators on each of the loops, that of the active power and 
that of the reactive power. Figure 5 shows the principle of the proposed method. 


> TE 
Psref.Qsre Linearized 


model 
of DFIG 
(DFIG+FBL) 





Eq(19) 





Figure 5. DFIG control scheme using proposed fuzzy-feedback linearization control (PFBL) 


The two most significant quantities to analyze the behavior of the system namely the power error 
ep Qand its variation dep, q, are chosen as two inputs of the regulator by Fuzzy logic controller FLC. 


{ ep, 7 Psret — Ps (45) 
ego, = Qsref g Qs 


Consider for each variable of measurement (the error and the variation of the error) 3 membership 
functions noted {N, EZ, P} with: Negative, About Zero, Positive presented in Figure 6. 


u(e), (de), u(u) 
N Z P 





-1 0 1 edeu 


Figure 6. Membership functions 


The rule bases of power controller’s are illustrated by the following Table 1: 


Table 1. Basis of fuzzy control rules 
e 








U N EZ P 

N N N EZ 
de EZ N EZ P 
P EZ P P 
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For the defuzzification of the output variables of the system, we use the method of the center of 
gravity (COG) [22]. 


_ SM nce) 
USP uGp (46) 


The proposed control is represented in Figure 7. 
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Figure 7. Fuzzy-feedback linearization control combined with high gain observer 


6. SIMULATIONS AND RESULTS 

In what follows we will study the performances (Reference tracking, disturbance sensitivity and 
robustness) for different controllers mentioned above conventional feedback linearization control (CFBL) 
and proposed Fuzzy-Feedback Linearization control (PFBL) using the rotor flux observer. 


6.1. Reference tracking 

This test consists in making a change in the active and reactive power setp values while maintaining 
the drive speed of the generator constant. Figure 8 illustrates the behavior of active and reactive powers 
stator. By examining this figure, it can be seen a good reference tracking with a less ripples for the proposed 
method compared with the conventional method. The decoupling between the two powers is ideally noted for 
the proposed method contrary to the conventional method we notice the appearance of a static error at times 
of step change (see Figure 9). In Figure 10, the observer's performance is illustrated when the estimated rotor 
fluxes and real fluxes are substantially identical for the dq axis. 
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Figure 8. Responses of active and reactive powers 
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Figure 9. Zoom responses of active and reactive powers 
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Figure 11. Zoom direct and quadrature flux. 


6.2. Disturbance sensitivity 


This test use to verify to what extent the measured powers stay at their reference when the grid 


voltage drops by 20% between 0.8s and 1.3s. The effect of this fault on the active and reactive powers of the 
machine is illustrated in Figure 12. 


The power references are correctly tracked, except for the presence of oscillations which are greatly 


increased during the fault, the power measurements show a significant deviation from the value of the 
reference and a reset time greater than 400 ms in the case of a conventional controller compared with the 
proposed controller which more effectively rejects the voltage drop effects (see Figure 13). The performances 
of the observer are verified by the simulation results presented in Figure 14. It can be clearly seen that this 
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voltage drop does not affect the observer because the estimated flux converge to their real values during the 


fault before stabilizing at steady state (see Figure 15). 
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Figure 12. Responses of active and reactive powers during stator voltage drop. 
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Figure 15. Zoom direct and quadrature flux during stator voltage drop 
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6.3. Robustness test 

Figure 16 shows the evolution of the active and reactive powers during a parametric variation of the 
generator when the resistance values are increased by 200% and the inductance values decreased by 60%. 
These results show that the excessive parametric variation caused a clearly degradation of the active and 
reactive powers curves with the appearance of a static error in the case of the conventional controller. The 
response of the system with the proposed controller remains insensitive to these variations. 
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Figure 16. Dynamic responses of active and reactive power under parametric variation of DFIG. 


7. CONCLUSION 

This work presents a feedback linearization control improvement for a DFIG-based wind system. In 
this context, we applied the fuzzy logic controller to the feedback linearization strategy. This is intended to 
improve the performance of the system particularly against parametric uncertainties and disturbance 
sensitivity that have an effect on conventional feedback linearization control used PI. From the comparative 
results between the two approaches presented, it can be concluded that intelligent controller have significant 
improvements to dynamic performances over conventional controllers. In addition, a part was reserved for 
the synthesis of a high-gain observer to reconstruct non-measurable rotor flux components for technical and 
economic constraints. This observer is tested by numerical simulation in combination with the adopted 
control that meets the assigned objectives. 











APPENDIX 
Table.2 Wind Turbine System Parameters 
Parameters value unit 
Nominal Power 1.5 MW 
Turbine radius 35.25 m 
Gearbox gain 90 
Stator Voltage 398/690 Vv 
Stator frequency 50 Hz 
Number of pairs poles 2 
Nominal speed 150 Rad/sec 
Stator resistance 0.012 Q 
Rotor resistance 0.021 Q 
Stator inductance 0.0137 H 
Rotor inductance 0.0136 H 
Mutual inductance 0.0135 H 
Inertia 1000 Kg.m2 
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